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INTRODUCTION 


One of the most common ways to characterize sound absorption and sound 
transmission in structures is through the use of sound power techniques (i,e. 
energy methods) in reverberation chambers. Sound power methods allow one to 
reduce an acoustic emission, absorption, or transmission problem to an 
equivalent electrical circuit problem with a power source and a network of 
complex impedances (impedances with resistive and reactive components). 

In order to use power flow methods in a reverberation chamber, a diffuse 
sound field is a necessary condition. Schroeder^ defines diffusion of a sound 
field at a point as "the angular distribution of sound energy flux in the 
plane wave expansion of the sound field at that point. If the distribution 
over the solid angle is uniform, we shall call the sound field at this point 
completely diffuse ." Rewording this, one might say that complete diffusion of 
the sound field at a point occurs when the time averaged acoustic intensity 
vector at that point is uniform in all directions. 

A substantial amount of research has been devoted to the study of sound- 
field diffusivity in reverberation rooms [see refs. 1-9]. Many different 
approaches to the problem have been investigated. The earlier studies [see 
refs. 2, 3] were primarily analytical attempts to describe room character 
since measurement equipment was primitive and computers were unavailable. 
Bolt^»^ took a statistical approach to the problem and tried to link the 
sound-field diffusivity to the "Frequency Irregularity" of the natural fre- 
quencies of the room. Later papers by Schroeder^ »^ , Doak** and others showed 
that the frequency irregularity is not related to the sound diffusivity 





in large rooms. The experimental portion of Schroeder's papers^ 9 ^ concen- 
trated on methods for measuring sound-field diffusivity at the wall surfaces 
of a room, A paper by Sepmeyer*’ used numerical schemes to show that there 
exists no truly "optimum" dimensions for a rectangular reverberation room, 
Doak^ provides a summary of the various theories dealing with the topic of 
room character. Some of the important conclusions reached in the various 
papers^’"^ are: 

(1) The -randomness of the acoustic field in a rectangular room with , 
smooth hard walls peaks at a particular frequency, and then decreases 
with increasing frequency [see refs. 3, 4, 7], 

(2) Above a certain critical frequency the statistical parameters of fre- 
quency response curves for all rooms are either identical or depend at 
most on reverberation time [ see refs, 4, S]„ 

(3) Averaging sound pressure over a frequency band containing a 
sufficient number of eigen frequencies appears to be equivalent to 
averaging over source and receiver positions as well [see ref. 4]. 

(4) There exist no truly "optimum" room dimensions for a reverberation 
room [see ref. 6]. 

The purpose of the research for the present paper was to determine the 
acoustic properties of the hard-walled acoustic facilities at NASA Langley 
Research Center to insure that the facilities meet the necessary requirements 
for the implementation of power flow techniques, 

THEORY AND INSTRUMENTATION FOR THE SOUND-FIELD DIFFUSIVITY TESTS 
A simple test for sound field diffusion was needed to test the quality of 
the sound-field of the various hard-walled acoustic facilities at NASA's 

LaRC. These acoustic facilities have irregular geometries making it difficult 
to calculate the eigen functions and natural frequencies of the rooms. 
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Consequently, an analytical study for determining room character was 
impractical. Although the literature^"^ provided useful information on the 
character of reverberation chambers, only one paper® presented a simple test 
for determining the sound-field diffusion in a room. 

Cook, et al® presents a simple experimental technique using two 
microphones for determining the sound-field diffusion in a room regardless of 
the room geometry. In this paper® it is shown that in a three-dimensional 
acoustically diffuse field, the cross-correlation coefficient between two 
microphones separated by a distance r is given by 

Pi2(0) = sin(kr)/kr , (7) 

where Pi 2 ( 0 ) is the cross-correlation coefficient (evaluated with zero time 
delay) and k is the wave number given by 


k = 


2irf 

Co 


( 8 ) 


where f is frequency and Cq is the speed of sound. Cook® also shows that 
for a two-dimensional diffuse field in which the acoustic intensity vector is 
confined to two dimensions, the cross-correlation coefficient between two 
microphones separated by a distance r is given by 

Pi2(0) = Jo(kr) (9) 

where Pi2(0) and k are defined by equation (7) and Jq is the 0^*^ order 
Bessel function of the first kind. The theory can be extended to show that in 
a one-dimensional acoustic field (e.g. the field in an infinitely long tube) 
the cross-correlation coefficient between two microphones spaced at a distance 
r will be 

Pi 2 ( 0 ) = cos(kr) . (10) 
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The cross-correlation coefficients are defined by the equation [see ref. 11] 


P12(t) 


J Rii(T)R22('r) 


( 11 ) 


where the subscripts 1 and 2 denote the two microphones and x represents the time 
delay. Mathematically the auto- and cross-correlation functions are given 
respectively by 


Rxx('t) = / x(t+x) dt and, (12) 

.00 

00 

Ryv^"^) " / X(t+x) y(t) dt . (13) 

-00 


The cross correlation coefficient, Px2(0) can be calculated using the auto- and 
cross-power spectral density functions [see ref. 11]. The auto- and 
cross-spectral density functions are given respectively by 


Gxx(f) = 2T[X(f) 

* X*(f)] and. 

(14) 

Gxy(f) = 2T[X(f) 

• Y*(f)] , 

(15) 


where T is the length of the time record and X(f), Y(f) are the complex finite 
Fourier transforms of the two signals given by 


X[f) 

T 


T 

/ 


0 


x(t) e^“^ dt 


(16) 


4 


Y(f) = - / y(t) dt 
T 0 


(17) 



where the asterisk (*) denotes the complex conjugate. With this information the 
cross-correlation coefficient between microphones 1 and 2 may be calculated using 
the equation 

I ReCGi2(f)] df 

Pi2(0) = ■ '■■■■■ -- (18) 

/ / Gii(f) df . / 622(f) df 

where Re[»] denotes the real part of a complex function and the limits of integra- 
tion depend on the maximum frequency of interest. 

For the experimental data in this paper, this method for calculating Pl2(0) 
was used. Since it is desirable to calculate pi2(0) as a function of kr (i.e. as a 
function of frequency f) equation (18) must be modified slightly to 


Pi2(T=0;f) 


Re[Gi2(f)] 

/ Gii(f) • G22(f) 


(19) 


The equipment used for the tests is shown in the block diagram of figure 1. 

Two microphones were placed in the room under test at a fixed distance apart, 

.3048 m (1 ft.), with their diaphrams in the same plane. Independently generated 
white noise signals drove each of the 2 (or sometimes 4) speakers in the room. 
One-third octave band equalizers were used to condition the input signals to the 
speakers to ensure that the noise field in the room was uniform broadband noise. 

The speakers were placed close to the corners of the room in an effort to excite as 
many of the room modes as possible. A dual channel Fast Fourier analyzer obtained 
the cross- and auto-spectral densities from the microphone signals over a 0-2000 Hz 
frequency range with a 5 Hz bandwidth. A desktop calculator computed pi2(0) as a 
function of frequency and plotted the results along with the theoretical curves. 

The experiment was performed with three microphone orientations. In each 
case, the microphone diaphrams were parallel to one of the three orthogonal 
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coordinate axis directions defined by the intersecting floor and walls of the 
room. The Fast Fourier analyzer took 400 ensemble averages for each 
measurement. This ensured that the random portion of the sampling error was 
no more than +; 5 percent (see ref. 11). 

Since instrumentation phase mis-match can detract from the accuracy of the 
measurements, the phase shift properties of the microphone-amplifier systems 
were carefully investigated prior to the experiments. The relative 
phase-match of the two channels due to the instrumentation is shown in figure 
2. The phase-shift is seen to be no more than 4° at 5000 Hz. 

REVERBERATION CHAMBER 

A series of sound-field diffusivity tests using cross-correlation 
coefficients were performed on the ANRL reverberation chamber (Bldg. 1208). 
Total volume and surface area of the chamber are 220 m^ (7769 ft^) and 229 m^ 
(2465 ft^) respectively. The floor of the chamber is smooth concrete and the 
ceiling and walls are segmented sections (splayed walls) of concrete. 

Prior to the sound-field diffusivity tests, two other properties of 
chamber were investigated. Assuming the chamber is a rectangular hard-walled 
enclosure, the first few eigen frequencies and the modal density as a function 
of frequency were calculated. The results of these calculations are given in 
table 1 and figure 3. The reverberation times of the chamber were measured in 
an earlier study^°. The results of the measured reverberation times of the 
chamber are shown in figure 4. 

Diffusivity tests on the reverberation chamber were performed using the 
equipment shown in figure 1. Four independently driven speakers positioned in 
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each of the four corners of the chamber were used to produce the uniform 
broadband noise field in the chamber. The microphones were placed 1.5 meters 
above the floor of the chamber in approximately the center of the working 
space of the room. The results of the cross-correlation measurements for the 
reverberation chamber are shown in figures 5 and 6. The narrowband (5 Hz 
bandwidth) cross-correlation coefficient measured over a 0-2000 Hz frequency 
range (0 to 12 kr range) is plotted in figure 5 for one of three different 
microphone orientations. The theoretical curves in figure 5 are the three 
dotted lines (calculated using equations (7), (9), and (10)) and the solid 
curve is the measurement data. It is evident from this figure that the mean 
value of the measurement data regresses to the theoretical curve for a diffuse 
acoustic field in three dimensions (sin[kr]/kr) . The standard deviation of 
the data about the theoretical curve, a, is .1366. Morrow® postulates that in 
order to obtain good agreement with Cook's theory it is necessary to average 
the data over a frequency band containing many room modes. If this condition 
is not satisfied large variations from the theoretical curves may be 
expected. A 50 Hz band moving average was performed on the data of figure 5. 
The result of this moving average is plotted in figure 6 along with the 
theoretical curve for a diffuse acoustic field in three dimensions 
(sin[kr]/kr) . Much better agreement with the theory is obtained by using this 
band average technique (a=.0649). Figure 6 indicates that the acoustic field 
in the reverberation chamber is acceptable for sound power measurements over 
the 100-2000 Hz frequency range. Similar results were obtained for the 
measurements with the other two microphone orientations. 
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TRANSMISSION LOSS FACILITY 


The ANRL transmission loss facility is a two-room acoustic test facility 
designed for noise reduction (NR) and transmissin loss (TL) measurements on 
light aircraft panels. It consists of two chambers, a source room and a 
receiving room, with an adjoining wall. The source room is used to create an 
incident sound field on the aircraft panel under test. The second chamber, 
called the receiving room, is used to measure the noise transmitted through 
the aircraft panel. The aircraft panel itself is positioned in a 1.23 m by 
1.525 m section in the adjoining wall between the two rooms. 

The source room of this facility measures 2.7 m x 2.9 m x 3.9 m (9‘ x 
9.5' x 12.75'). Total volume and surface area of the room are 30.9 m^ (1090 
ft^) and 59.7 m^ (645 ft^) respecti vely. The room has a tile floor, plaster 
board walls and a suspended acoustic tile ceiling. Assuming the room is a 
rectangular hardwalled enclosure, the first few eigen frequencies and the 
modal density as a function of frequency were calculated. These calculations 
are given in table 2 and figure 7. Measured and calculated octave band 
reverberation times are shown in figure 8. 

Reverberation time was calculated using the Mi 1 1 ington-Sette^^ equation 
given by 




.049V 

I S.ln(l-a.) 


( 20 ) 


where Teo is reverberation time, V(ft^) is the room volume, Si(ft^) is 
surface area of the i^^*^ surface in the room, and a-j is the acoustic random 
incidence absorption coefficient of the ii^*^ material in the room. The 
values of V and Si were measured values and the ai values were obtained 
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from reference 13 and measurements. The values of a-j used in the calcula- 
tions are given in table 3= 

Measurements of reverberation time were taken several times for each 
octave band at a single location in the room using a commercially available 
reverberation time meter. The results of the measurements were averaged into 
a single reverberation time for each octave band and the results plotted in 
figure 8. This figure shows that the measured and predicted values of the 
reverberation time agree in the lower frequency bands only when the airspace 
above the suspended ceiling is taken into consideration in the calculations. 
This indicates that below 1 kHz, the airspace above the ceiling has a 
significant effect on the reverberation time of the room. 

Diffusivity tests on the souce room were performed using the equipment 
shown in figure 1. Two independently driven speakers, positioned in corners 
of the room, opposite the aircraft panel, were used to produce the uniform 
broadband noise field in the chamber. The microphones were placed 1.5 meters 
above the floor of the room in approximately the center of the working space 
of the room. An all aluminum aircraft panel with stringers and frames separ- 
ated the source room and receiving room for the duration of the diffusivity 
tests. An example of the results of the cross correlation measurements are 
shown in figures 9 and 10. The narrowband (5 Hz bandwidth) cross correlation 
coefficient measured over a 0-2000 Hz frequency range is plotted against the 
theoretical curves in figure 9 for one of three different microphone 
orientations. The measurement data of figure 9 exhibit a significant amount 
of scatter about the theoretical curve for a three dimensional diffuse field 
(sin[kr]/kr) . The standard deviation of the data about the theoretical curve, 
a, is .2954. A 50 Hz moving band average of the measurement data was 
performed and is shown in figure 10, This figure shows that the data deviates 
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considerably from the theoretical curve (a=.1529). Comparison of figure 10 
(source room) with figure 6 (reverberation chamber) indicates that the quality 
of the sound field in the reverberation chamber is superior to that of the 
source room of the transmission loss facility as expected. 

In an effort to improve the quality of diffuseness in the source room, 
the suspended acoustical tile ceiling was removed. This increased the height 
of the room by 1 meter so that the total volume and surface area of the room 
are 41.9 m^ and 72.9 m^ respectively. The ceiling of the room, without the 
acoustic tile, consists of irregularly shaped steel reinforced concrete. 

Reverberation time was measured, as before, with the source room in its 
new configuration. These data were compared with the previously measured 
reverberation times. The comparison is shown in figure 11. This figure shows 
that the reverberation time of the source room was increased substantially as 
expected by removing the acoustical tile ceiling. 

Cross-correlation diffusivity tests were performed on the source room 
with its new ceiling configuration using the same techniques discussed pre- 
viously. The results are shown in figures 12 and 13. The standard deviation 
of the narrowband measurement data about the theoretical curve (sin[kr]/kr.) in 
figure 12 is o=.2969. Comparison of figures 10 and 13 indicates that the dif- 
fuseness of the source room was improved substantially by removing the acous- 
tical tile ceiling (o=.1277). The enhanced quality of the sound field in the 
room is particularly evident in the frequency range of 500-1500 Hz (3 to 9 kr). 

The receiving room of the ANRL transmission loss facility measures 4.06 m 
X 2.9 m X 2.74 m (13-1/3' x 9-1/2' x 9'). Total volume and surface area of 
the room are 32.3 m^ (1140 ft^) and 61.7 m^ (664 ft^) respectively. The room 
has a tile floor, plasterboard walls and a rigid-backed acoustic tile ceiling. 
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The receiving room is dynamically isolated from the building which houses the 
facility and has double concrete walls with a double door. Measurements of 
flanking transmission loss (for sound entering the receiving room via paths 
other than the panel) indicate that the receiving room may be used to measure 
sound levels that are 80 dB lower than the sound levels in the source room. 

The ambient noise level of the receiving space is less than 20 dB which is gen- 
erally below the noise floor of the condenser microphones used for measure- 
ments. 

Assuming that the receiving room is a rectangular hardwalled enclosure, 
the first few eigen-frequencies and the modal density as a function of fre- 
quency were calculated. These calculations are given in table 4 and figure 14. 
Measured and calculated octave band reverberation times are shown in figure 15. 

For many applications, such as the measurement of transmission loss, it 
may be desired that the receiving room have a very high absorptivity. If trans- 
mission loss is being measured by an array of microphones in the receiving space 
or an acoustic intensity meter (see ref. 14), it is necessary to increase the 
acoustic absorption in the receiving room to a level which approaches the absorp- 
tion of an anechoic space. An attempt was made to increase the absorption of the 
receiving room by placing six 1.2 m x 2.5 m x .10 m thick fiberglass panels in 
the room. The fiberglass panels were placed parallel to the wall adjoining the 
source and receiving rooms and were evenly spaced 60 cm apart. Reverberation 
times in octave bands were measured, as before, with the receiving room in this 
new configuration. A comparison of the results is shown in figure 16. This 
figure demonstrates that the absorption in the receiving space can be increased 
dramatically by the inclusion of relatively few fiberglass panels. 
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CONCLUSIONS 


The results of this indicate that it is possible to create a diffuse sound 
field in the ANRL reverberation chamber over a 100-2000 Hz frequency range. 
These results were for a specific system (viz, four independently driven 
speakers in each corner of the chamber) and data were taken only near the 
center of the room. It is possible that the quality of the sound field in this 
frequency range as well as higher frequencies will differ markedly depending on 
the system used to drive the sound field and on the location of the source and 
receiver in the chamber. 

The quality of the sound field in the source room of the transmission loss 
facility is acceptable for measurement of transmission loss over a 200-2000 Hz 
frequency range when using two independently driven speakers to produce the 
incident sound field. Random incidence transmission loss measurements in the 
frquency range below 200 Hz would be suspect due to the low modal density of 
the source room in this frequency range. 

The absorption of the receiving room of the transmission loss facility can 
be increased considerably by incorporating fiberglass panels. 
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TABLE 1 


REVERBERATION CHAMBER 
CALCULATED EIGEN FREQUENCIES 


MODE NUMBERS 

M M M 

' 'x y z 

1 0 0 

0 1 0 

1 1 0 

0 0 1 

2 0 0 

1 0 1 

0 1 1 

2 1 0 

1 1 1 


FREQUENCY, Hz 

21.6 

28.1 

35.5 

40.2 

43.3 
45.7 
49.1 

51.6 

53.6 


0 


2 


0 


56.3 



TABLE 2 


TRANSMISSION LOSS FACILITY - SOURCE ROOM 
CALCULATED EIGEN FREQUENCIES 


MODE NUMBERS 
^ M 

1 0 0 

0 1 0 

0 0 1 

1 1 0 

1 0 1 

0 1 1 

2 0 0 

1 1 1 

2 1 0 


z 


FREQUENCY, Hz 


44.13 

59.23 

62.52 
73.86 

76.52 
86.12 
88.26 
96.77 

106.29 


2 


0 


1 


108.16 



TABLE 3 


ABSORPTION COEFFICIENTS, g 


MATERIAL 

OCTAVE BAND CENTER FREQUENCY, 

Hz 


125 

250 

500 

1000 

2000 

4000 

PLASTER WALL 

.013 

.015 

.02 

.03 

.04 

.05 

♦ACOUSTIC CEILING 
TILE 

.1 

.15 

.25 

.3 

.4 

.4 

tACOUSTIC CEILING 
TILE W/ AIRSPACE 
ABOVE CEILING 

.3 

.45 

.50 

.4 

.4 

.4 

GLASS 

.35 

.25 

.18 

.12 

.07 

.04 

CONCRETE 

.01 

.01 

.01 

.02 

.02 

.02 


* Values for the acoustic ceiling were measured using the 
NASA Langley Research Center's impedance tube facility. 

t Corrections for the airspace above the ceiling were 
taken from Fig 10.13 of reference [12]. 






TABLE 4 


TRANSMISSION LOSS FACILITY - RECEIVING ROOM 
CALCULATED EIGEN FREQUENCIES 


M 

1 

0 

1 

0 

2 

1 

0 

2 

1 

0 


MODE NUMBERS 




0 

1 

1 

0 

0 

0 

1 

1 

1 

2 


0 

0 

0 

1 

0 

1 

1 

0 

1 

0 


FREQUENCY. Hz 


42.21 

59.23 

62.52 

72.73 

75.43 

84.42 

86.12 

95.91 

103.13 

105.15 



ROOM UNDER TEST 



Figure 1.- Instrumentation block diagram for the diffusivity tests. 
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Figure 2.- Instrumentation phase mis-match for the measurement system. 
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Figiire 3.- Calculated modal density of the reverberation chamber. 
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Figure 4.- Measured reverberation time in the reverberation chamber. 
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Figure 6,- Band averaged cross correlation coefficient measurements of 

diffusivity in the reverberation chamber. 
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Figure 7.- Calculated modal density of the source room. 
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Figure 8.- Reverberation time in the source room. 
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Figure 9.- Results of the cross correlation coefficient measurements of 

diffusivity in the source room. 
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Figure 10.- Band averaged cross correlation coefficient measurements of 

diffusivity in the source room. 
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Figure 11.- Increase in reverberation time in the source room 
with the ceiling modifications. 
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Figure 14.- Calculated modal density of the receiving room. 



RECEIVING ROOM 


y 


4 


MEASURED 



160 250 400 630 1k 1.6k 2.5k 4k 6.3k 

Frequency (Hz) 



Figure 15.- Reverberation time in the receiving room. 
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Figure 16.- Decrease in reverberation time in the receiving room 
with the addition of fiberglass panels. 
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